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Abstract 

Continuous unitary transformations can be used to diagonalize or approximately diagonalize a given 
Hamiltonian. In the last four years, this method has been applied to a variety of models of condensed 
matter physics and field theory. With a new generator for the continuous unitary transformation proposed 
in this paper one can avoid some of the problems of former applications. General properties of the new 
generator are derived. It turns out that the new generator is especially useful for Hamiltonians with a 
banded structure. Two examples, the Lipkin model, and the spin-boson model are discussed in detail. 
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1 Introduction 



The diagonalization of a given Hamiltonian is one of the important goals in any quantum mechanical 
problem. Apart from few explicitly solved models it is only possible approximatively or numerically. 
Four years ago, Wegner §\§ proposed flow equations for Hamiltonians to bring a given Hamiltonian closer 
to diagonalization. Technically, the approach is based on a continuous sequence of infinitesimal unitary 
transformations applied to the Hamiltonian. The infinitesimal unitary transformations are chosen so that 
the off-diagonal matrix elements become smaller. In principle the flow equations proposed by Wegner 
finally yield a diagonal matrix, except for the situation close to resonances, where few off-diagonal matrix 
elements may still be large. This will be made more precise later. Unfortunately it turned out that for a 
realistic Hamiltonian in an infinite-dimensional Hilbert space, the situation may be more difficult. Wegner 
discussed as an example interacting electrons in one dimension. The first problem in such a system is, that 
during the flow higher interactions are generated. Wegner solved this problem by passing to an ^-orbital 
model, where, in the large rc-limit, the flow equations can be closed. The second problem is, that even 
then it is not possible to diagonalize the Hamiltonian completely. Instead he used a block-diagonalization. 
With this modified approach he was able to solve the flow equations. 

Later Wegner's flow equations have been applied to various models of condensed matter physics. 
A class of models where the flow equations yield very accurate results are dissipative quantum systems 
[§, [|]. It was also possible to treat the electron-phonon problem with this method [Q, ^J|, for which one can 
obtain accurate values for the transition temperature The reason why the approach is so successful 
is that it provides a consistent renormalization scheme for Hamiltonians. This has first been pointed out 
by Glazek and Wilson [[7|], who developed an essentially equivalent method a few months later. Their 
method has been used to treat problems from quantum chromo dynamics H |9p. 

Flow equations are a useful tool to treat systems with various energy scales and to renormalize a 
given Hamiltonian, but it is difficult to apply them to a finite matrix or to a Hamiltonian with a simple 
structure. Although it may be possible to choose the transformation in such a way that some of the higher 
interactions are not generated (this has been done in the treatment of dissipative quantum systems in 
[§]), one cannot avoid in general that the Hamiltonian looses its initially simple structure. This point 
has been discussed in detail by Richter [10]. He applied the flow equations to a simple model, the spin- 
boson model, and used it to test various truncation schemes. The aim of the present work is to propose a 
different set of flow equations that has the property that a band diagonal matrix or Hamiltonian keeps its 
band diagonal structure during the flow. 

In the following section I introduce the new generator for the continuous unitary transformation and 
I show that it can be used to diagonalize a given Hamiltonian. The main property of the new generator 
is that band diagonal Hamiltonians remain band diagonal. I derive some properties of the new flow 
equations. 

In section 3 and 4 I apply the flow equations to two simple models: The Lipkin model and the spin- 
boson model. The Lipkin model has a finite Hilbert space, the Hamiltonian can be written in the form of 
two tridiagonal matrices. The spin-boson model can as well be written as two tridiagonal matrices, but 
the Hilbert space is infinite. The aim of this paper is to show that the new set of flow equations is useful 
if one wants to deal with simple finite or infinite matrices. Therefore I only derive some results for the 
spectra of these models, but I do not discuss the physics of these models in detail. In both cases the reader 
may consult the references for the physical background of the models. For the Lipkin model, all results 
presented here are well known and have been obtained by various other methods. For the spin-boson 
model I derive a formula for higher eigenvalues that has been derived so far only within a first order 
perturbational treatment. The new result is that this formula has a much wider range of validity. 

Section 5 contains the conclusions together with a critical discussion of the possible range of appli- 
cability of the new flow equations. 
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2 Generalities 



In this section I deal with a Hamiltonian that is given by a finite or infinite, real, symmetric matrix 

H — (hntn)i hnm — h mn . (2-1) 

An extension to complex, hermitian or normal matrices is easily done. In general, flow equations for a 
Hamiltonian are constructed by a continuous unitary transformation written in a differential form, 

dH . 

^ = m,H] (2.2) 

r| is the generator of the infinitesimal unitary transformation, it is an anti-hermitian operator that depends 
on H and therefore implicitly on the flow parameter I. Wegner proposed to choose r\ = [Hd,H], where 
is the diagonal part of the Hamiltonian. With this choice of r\ one can show that r\ —> for £ — > oo. 
The flow equations yield a final matrix with the property that h nm (h nn — h mm ) = 0. This means that either 
the off-diagonal matrix elements h nm vanish or that for a finite h nm the difference of the corresponding 
diagonal matrix elements h nn — h mm must vanish. Thus the Hamiltonian is diagonalized except for some 
possible resonances. As already mentioned in the introduction this choice of r\ has one disadvantage: If 
the initial matrix has a simple structure, it looses this structure for finite t. This becomes clear if one takes 
a band diagonal Hamiltonian, i.e. 

h nm (i = 0) = if \n-m\ >M (2.3) 
Taking 

Tl = (Tbim) (2-4) 
one obtains flow equations for the matrix elements 

dh nm 



di k 



YS^nkhkm - h nk r[hn) (2.5) 



which in general do not conserve (2.3). Therefore I propose the new generator 

f\nm = -Tlmn = sign(n - m)h nm , t\ m = (2.6) 

With this choice the flow equations for the off-diagonal matrix elements are 

dh 

—jj- = -sign(n - m) (h m - h mm )h nm + £ (sign(« -k) + sign(m - k))h nk h km . (2.7) 

Due to the sum of the two sign-functions the second term vanishes if |n — m\ > M. This shows that the 
new generator preserves the band diagonal structure. For the diagonal matrix elements one obtains 

-^p = 2 £ sign(« - k)h nk h kn (2.8) 

The main question is now whether the new choice of r\ can be used to diagonalize the Hamiltonian. This 
is indeed the case for finite matrices (1 < n < N) or semi-infinite matrices (1 < n). For the sum of the 
first r diagonal matrix elements one obtains the differential equation 

d r r 

-77 52 h nn = ~ 2 YjYj h "k h kn < (2.9) 
n=lk>r 
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This quantity decays as a function of £. I assume that H is bounded from below. Then Y!n=\ hnn is bounded 
from below by the sum of lowest r eigenvalues of H. Therefore its derivative must vanish in the limit 
£ —y oo, i.e. 

]imh n khkn = (2.10) 

Furthermore, if h nm tends to zero, one must have sign(« — m) (h nn — h mm ) < for sufficiently large values 
of I. Thus the diagonal matrix elements are ordered for large I. 



Some properties of the flow equations with the new choice (2.6) of the generator are: 

• The final matrix is diagonal, even if it contains degeneracies. 

• The matrix remains banded, h nm = if \n — m\ > M. 

• h nn {oo) > h mm (°°) for n > m and for irreducible matrices. If the matrix is reducible, each of the 
irreducible blocks can be treated separately. 

• The asymptotic behaviour of h nm for large I is known: h nm oc exp(— \h nn {o°) — h mm (°°)\£). 

The last property follows directly from the flow equations. For large values of £, the diagonal matrix 
elements are in the correct order and no level crossings occur for larger values of £. The second term in 



( p77[ ) falls of faster than oc exp(— \ h nn (oo) — h mm (°o)\£) so that the asymptotic behaviour is determined by 
the first term. For Wegner's choice of T|, one cannot exclude non-vanishing off-diagonal matrix elements 
due to degeneracies of diagonal matrix elements, and the asymptotic behaviour of the matrix elements for 
large £ is less clear. 

Let me now discuss the application of the new flow equations to two examples. 

3 Example 1: Lipkin model 



The Lipkin model [jll|] is a toy model of nuclear physics that describes in its simplest version two shells 
for the nucleons and an interaction between nucleons in different shells. It serves as a standard example 



for testing of various approximations. Recently Pirner and Friman Q12Q applied flow equations to this 
model. As usual new interactions are generated and they used a suitable truncation to close the flow 
equations. They showed that for a large number of particles N the truncated flow equations yield the 
exact result whereas for small N deviations occur. Furthermore the flow equations are only applicable for 
small couplings. For larger couplings the model shows a transition from a state with the same symmetry 



properties as the Hartree-Fock state to a state with different properties Q13[. I will come back to this point 
later. 



A suitable representation for the Lipkin model is in terms of pseudo-spin operators [|12|]. 

H(£ = 0) = ^ J z + V (J 2 + +J 2 ). (3.1) 

H commutes with J 2 = j\ + i(7 + /_ +/_/+) and J 2 has the eigenvalues 7(7+ 1) as usual. The pseudo- 
spin operators form the usual angular momentum algebra 

[7 Z ,7±]=±7 ± , [7+,7_]=27 z . (3.2) 

It can easily be shown that in the basis where J 2 and J z are diagonal the Hamiltonian decays into two 
tridiagonal matrices. 

H = (h nm ) n ^ m= Q j or j-i or 7-1/2 (3-3) 
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The dimension of the matrices depends on /. If 2/ is even, one of the two matrices has dimension /, the 
other 7+ 1, if 2/ is odd, both matrices have the dimension / + 1/2. The matrix elements are 

h nn = e„, h nn+ \ = S„, h nm = if \n-m\>l (3.4) 

with 

£„(0) = ^„(-/ + 2«) (3.5) 



S„(0) = V oy /J(J+ !)-(/- 2n)(J -In- \)yJJ(J+ 1) - (/ - In - 1)(/ - In - 2) (3.6) 



or 



e„(0)=£ (-/ + 2n + l) (3.7) 



8„(0) = VbV- / (- / + 1 )-(- / - 2 «- 1 )(- / - 2 «- 2 )v / - / (- / + 1 )-(- / - 2 "- 2 )(- / - 2 «- 3 )- ( 3 - 8 ) 
The flow equations are in both cases 

^ = -28£ + 28g_i (3-9) 

-^ = -8„(e n+1 -e„). (3.10) 

A first possibility is to solve these equations iteratively. One can start with the ansatz (£) = z n (0) and 
§1°^) = 8„(0)exp(— 2^o^)- Inserting these expressions on the right hand side of the flow equations yields 
a first iterative solution, which can again be used to obtain the next iterative solution and so on. When 
one uses this procedure it may be useful to write the flow equation for 8„ in the form = — £„ + i +£„. 
This procedure reproduces simply perturbation theory, which works well for small Vq and not too large / 
(i.e. not too large particle numbers). 

A simple non-perturbative solution can be obtained in the limit of large /; this corresponds to the 
limit of a large particle number. The two different cases for the initial conditions above yield 

8, 1 (0) 2 -8, 1 _ 1 (0) 2 = 32y () 2 / 2 (« + i)(l + O(l/7)), (3.11) 
or 

8„(0) 2 - S^O) 2 = 32V 2 / 2 ( n + 1)(1 + 0(1//)). (3.12) 
With the ansatz 

£„(£)= a(£)n + b(£) (3.13) 
the flow equations can be written as 

8 n (£) = f(£)8 n (0) (3.14) 

where 
df 

and 

^ = -64V 2 / 2 / 2 . (3.16) 
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For the last two equations the quantity 

a 2 -64V^J 2 f (3.17) 
is conserved. Since / — » for I — > oo, this yields directly 



«H = V 4 ^o-64V r 2 / 2 (3.18) 

This solution exists for 4/Vb < ^o- Taking the first case from above one has 

db[ 1 da 
"jF = 4dl 
with the solution 

1,. i 



(3.19) 



bi H = - (J+ -)§o + - ^ - 64V 2 / 2 (3.20) 



For the second case one obtains 

dbi 3 
"jF ~ 4d£ 
and 



(3.21) 



^(oo) = - {J+ ^ + y 4 z2_ 64v 2j2 (3 . 22) 

This yields directly the approximate spectrum 

1.1,, 1 



e„i, 2 = y^g - 64V 2 / 2 (n + - ± - ) - (/ + -)£o (3.23) 
and the gap between the ground state and the first excited state 



602-60! = J^-lGfiJ 2 . ( 3 - 24 ) 



This result has been obtained by Pirner and Friman Q12| ] as well. It is also well known from RPA 



Nevertheless the above formulation of the flow equations has certain advantages: 

1. It is quite easy to obtain perturbation theory using the flow equations. With the conventional for- 
mulation of the flow equations this in principle possible, but one has to introduce many higher 
interactions if one wants to obtain higher orders in perturbation theory. In the present formulation 
the flow equations (3^, 3.10| ) are closed and the iterative solution is easily constructed. 



2. The formulation of the flow equations is not restricted to 4JVq < ^o- Although the perturbative 
solution and the approximate solution for large / shown here are limited to this regime, the flow 



equations (3.9, 3.10) can be solved (at least numerically) for 4JVq > t,o as well. 



3. The flow equations in the present form may be used to derive a systematic 1 //-expansion. To do 
this one has to use a polynomial ansatz for £„ as a function of n instead of the linear ansatz above, 
and one has to take higher orders in 8„ — 8 2 _ j into account. 

Since the aim of the present paper is only to show that the new proposal for the continuous unitary trans- 
formation is useful if one wants to treat band-diagonal Hamiltonians, I do not follow the lines suggested 
in these points. 
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4 Example 2: Spin-boson model 



As a second example I discuss the spin-boson model described by the Hamiltonian 

H(£ = 0) = --<5 x + -c z (b + b } ) + (ab } b. (4.1) 

It has a wide range of possible applications, especially in atomic physics where the spin describes a 
two level atom that is coupled to a e.g. laser field. Due to its long history there exists an enormous 
amount of work that has already been published on this model, so that it is impossible to review or cite 
all these papers. A good overview may be found in the paper by Graham et al. [ff3|. More recently 



this model has been discussed in connection with quantum chaos [|15|, |16J]. Together with the usual flow 
equations the model has been used to test several approximation schemes [|K]]. It turned out that the 
ground state and the low lying excited states as well as dynamical properties can be calculated very 
accurately using traditional flow equations. If one is interested in quantum chaos, an accurate knowledge 
of high eigenvalues is necessary. 



The Hamiltonian (4.1) can be written as two tridiagonal infinite matrices. The flow equations are 



therefore the same as for the Lipkin model, 

^ = -2S£ + 2J5£_i, (4.2) 
and 

^ = -5 n (e n+ i-e„), (4.3) 
at 

but with different initial conditions: 



5„(0) = -Vn~+\, (4.4) 

6„(0)=n(0±(-l) n |. (4.5) 
It is very easy to solve these equations for A = 0. One obtains 

£„ = «co + £o, 8o = — — (1 — exp(— 2o>£)) ; 8„ = -y/n+ lexpf— (o£). (4.6) 
4co 2 

In principle it is possible to use this solution to obtain an expansion for small A. A perturbative treatment 



for small A has been given to first order by Graham et al. []14[], and the flow equations yield the same 
result. Therefore I will not reproduce this solution here. A second possibility is an iterative solution for 
small X. It yields a perturbative solution valid for small n. Similarly ordinary perturbation theory is valid 
only if n is small (n <C (co± A) 2 /X 2 ). Instead I try to obtain an asymptotic expression for £„ that is vaild 
for large n. To do this, I make the following ansatz. 

e„ = raffl- £-(1 -ex V (-2(o£)) ± (-l) n ^f n (£), (4.7) 
4(0 2 

and 

52 = ^(« + l)exp(-2(o^)±(-ir^). (4.8) 
This ansatz yields flow equations for /„ and g n 

^ = -2(g n + g n _ 1 ) (4.9) 
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^ = ^(„+l)exp(-2ro£)(/ n+1 +/„)-2 ffl& TA(-l)^(/ B+1 +/„) (4.10) 
It is now useful to introduce the new variable 

x = 1 - exp(-2co£) (4.11) 

instead of the flow parameter I. In the following I take /„ and g n as functions of x. The flow equations 
are rewritten as 

(0(1 -xy^- = -g„-g n -i, (4.12) 
ax 

2(0(1 -x)^ = y (n+ 1)(1 -jc)(/ b+ i +/„) -2a>g„ TA(-l)%(/„ +1 +/„), (4.13) 
with the initial conditions 

/»(0) = 1, *„(<>) = 0. (4.14) 

These equations are still exact. For large n one has g2 n +2 ~ g2« and g 2ra +i ~ gin-\- As a consequence one 
obtains /2 n +i ~ /2 n - This yields 

0,(1 _ X )^L « 0,(1 - x ) d Jf^ = - g2n - g2n+l (4.15) 
ax ax 

2co( 1 - x) ^ « 2X 2 « ( 1 - x)/ 2n - 2(o^ 2 „ t 2Ag 2n f 2n (4.16) 

2(0(1 -^)^^ ~ 2X 2 «(1 -x)/ 2 „ -2(0 52 „ +1 ±2A 52n+1 / 2n (4.17) 

Taking the derivative of the first equation, and using the sum of the second and the third to express the 
derivative of g 2n + g2n+\ by fin I obtain for large n (i.e. 2n + 1 « 2n) 



The general solution of this differential equation can be expressed using Bessel functions 

aJ x y/nil-xfj + bY x {^L^/ n {\-x)\ (4. 19) 



/„(*) = \/\-x 



Using the initial conditions f„(0) = 1 and f n (0) = 0, one obtains 
a = Ti—y/nYo — \/n 

b = -k-x/EJq — y/n) (4.20) 



(0 V (0 



This yields the complete solution for £„(£). I am interested in the limit i — ► °°, which corresponds to 
x = 1. Using the behaviour of the Bessel function for small arguments one obtains 



-/o(f^) 



/n(l) = -/o T7V« (4-21) 
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This yields the final expression for the eigenvalues for large n 
4co 



n(0 



Using the asymptotic behaviour of the Bessel function, one obtains 



e„(£=~) = Mo-^=F(-l)"^Aj 



COS 



2X 
CO 



(4.22) 



(4.23) 



Except for a misprint (the factor n 1//4 is missing in the last term) this expression coincides with the result 
in Graham et al. performed a first order perturbational treatment in A and expanded the result to 



obtain ( 1.23 ) for large n. Thus, in their approach, the validity of ( 4.23| ) is unclear. 

How accurate are ( |4.22| , [4.23 )? A first condition is obtained from a consistency check of the above 
assumption f n +i(x) ~ /«(•*)• This assumption must be true for all x and all X/(a. Due to 2\Jn + lX/co « 
2y/nk/(a + X/ (co-y/w) one must have 

X 



< 1 



(4.24) 



But this is not the only condition one needs. A similar consistency check has to be done for g n . It is more 
complicate since it depends on A. Fortunately there is a more simple possibility to determine the range 
of applicability for ( |4.22 , (4.23 ), which is equivalent to the consistency check proposed above. From the 
general considerations in section 2 one knows that the flow equations yield e n (£ = °°) < £„+i(^ = °°). 
This means that one must have 

' 2X ^ • f2X ^Tl) . (4.25) 



co > 



Jo 



CO 



Jo 



CO 



For large j-\Jn this yields the condition 



A 
2co 




X 1 



7IC0 w 3 / 4 



< 1 



(4.26) 



( 4.24 ) and (126) show that within a wide range of parameters ( 4.22| , |4.23 ) are applicable. But one should 
be careful using ( [4.23 ) for small X. (4.23) has a relative error O(^j^) compared to (4.22) so that for 
small X ( 4.22j ) yields better results. For X = 0, ( 4.22| ) yields the exact solution (which is of course trivial), 
whereas (4.23) is not defined. Only in the limit X — > and AX/co 2 = const it is not possible to apply 



(4.22). But this regime can be treated with perturbation theory for small X. 



The accuracy of (4.22) is very high, even for small n and large A. It can be tested numerically if one 
compares ( 4.22 ) with the exact numerical solution of the flow equations. This is done in figure 1. We 



show the relative error of the asymptotic formula ( 4.22 ) plotted as a function of A/co. The coupling is 
X/co = 4.0. For small n the error is relatively large (up to 2% for n = 10), whereas it is smaller for larger 
n. The parameters are chosen so that for n = 10 the left hand side of ( 4.24 ) equals 1.3 and is therefore too 
large. The left hand side of ( 4.26| ) equals 0.5 for the worst case shown, i.e. n = 10 and A = 5(0. A relative 



error less than 0.1% is obtained if the left hand side of (126) is less than 0.1 and for n^X 2 /® 2 . 

The main result of this section is clearly the expression for z n {£ = oo) given above. But one does not 
only know the eigenvalues of the Hamiltonian, the complete flow ( |4.19| ) is known as well. This allows to 
reconstruct the continuous unitary transformation that diagonalizes the Hamiltonian. Therefore it is even 
possible to calculate other observables in the transformed basis. In principle it is possible to combine 
this result with the results of Richter [[R]], which are very accurate for small n. Then one should be able 
to obtain accurate values for dynamical correlation functions of the model. But the goal of the present 
calculations was to show that the new flow equations proposed in this paper are indeed useful if one wants 
to deal with band-diagonal matrices or Hamiltonians. Therefore I do not follow this line. 
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0.020 




Figure 1 : The relative error of the asymptotic expression for the eigenvalue compared to the exact result 
as a function of A/co for X/iO = 4.0. The short dashed line is n = 10, the long dashed line n = 15 and the 
solid line n = 20. 



5 Conclusions 

The flow equations proposed in this paper seem to be useful if one wants to solve a given model. The 
two examples showed that solutions can be obtained at least approximatively or in some limits. Some 
of the advantages of the new flow equations have already been pointed out: 1. The matrix is finally in 
a diagonal form, even if degeneracies occur. 2. If the Hamiltonian has initially a banded structure, this 
form is preserved. 3. The asymptotic behaviour of the equations is known. Another advantage is that 
with the new flow equations one can obtain accurate results for higher eigenvalues as well. But there are 
also some disadvantages, which may be important in other systems: 

1. Although asymptotically off-diagonal matrix elements decay faster if the difference of the corre- 
sponding diagonal matrix elements is larger, the new generator does not separate different energy 
scales automatically. For small I it is not guaranteed that off-diagonal matrix elements corre- 
sponding to large energy differences decay fast. Therefore the new flow equations do not provide a 
renormalization scheme. 

2. An important property of a Hamiltonian one wants to solve using the new flow equations is that it 
must have a pure point spectrum. Hamiltonians with continuous spectra cannot be treated that way. 
It is for instance not possible to apply the new generator to a dissipative quantum system. 

3. If the diagonal matrix elements of the initial Hamiltonian are not in the correct order (i.e. h nn < h mm 
if n < m) the flow equations will reorder the diagonal matrix elements. This may cause a problem, 
because in such a case the analytical treatment of the flow equations becomes more difficult. 

The first two points mentioned above show that the new generator cannot be applied successfully to the 
kind of problems that have been treated so far using flow equations. In all these problems one has several 
different energy scales and continuous spectra, and one needs a renormalization scheme to obtain useful 
results. In this sense the new generator provides a complementary set of flow equations. It can be applied 
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to problems that cannot be treated with the original flow equations, but problems that can be treated with 
the original flow equations are not within the range of possible applications of the new scheme. 

What kind of problems can be treated with the new flow equations? It became already clear that 
Hamiltonians with a banded structure are good candidates. But even for Hamiltonians without such a 
structure the present approach may be useful. If one has for instance a problem for which all diagonal 
matrix elements are of the same order of magnitude, the usual flow equations are difficult to apply. The 
reason is that differences of diagonal matrix elements are small and that therefore the flow is very slow. 
This does not happen with the new flow equations, since the flow of the diagonal matrix elements is only 
determined by the magnitude of the off-diagonal matrix elements. A class of possible candidates are 
therefore disordered systems. 

In general one can say that the new flow equations can be applied to single or few particle systems. It 
may therefore be useful in nuclear or atomic physics. 
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